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ABSTRACT 


s 

In  r.'ost  approaches  to  modelling  the  rainfall/runoff 
process,  a  spatially  lumped  description  of  precipitation  has 
been  assumed  adequate  for  modelling  the  important  aspects  of 
catchment  response.  However,  theories  of  catchment  hydrology 
as  well  as  some  recent  modelling  studies  suggest  that  spatial 
variability  in  precipitation  may  be  important  in  determining 
the  characteristics  of  streamflow  hydrographs.  Data  frcm  two 
intensive  rainfall  recording  experiments  in  Illinois  have 
been  used  to  examine  the  effects  of  rainfall  pattern  on  stream 
hydrographs  for  summer  convective  storms.  A  threshold  analysis 
was  used  to  distinguish  storms  of  markedly  different  pattern. 

A  mixed  deterministic/stochastic  modelling  procedure  was  used 
to  determine  the  length  of  record  required  to  differentiate 
the  hydrograph  characteristics  resulting  from  storms  of 
different  patterns.  It  was  found  that  differences  in  peak 
timing  were  highly  significant  but  that  differences  in  the 
distributions  of  peak  flow  and  stormflow  volumes  were 
generallv  insianif icant  even  given  a  long  period  of  record. 

/V 


I  NTRODUCTION 


accurate  assessment  of  the  spatial  pattern  of 
precipitation  delivered  to  the  surface  of  a  catchment  may  be 
important  for  rainfall /runoff  modelling  from  two  viewpoints. 
First,  derivation  of  a  useable  estimate  of  the  total  volume 
of  storm  precipitation  for  a  basin  may  require  that  the 
spatial  distribution  of  rainfall  be  taken  into  account.  The 
question  of  spatial  variation  in  this  case  is  related  to  the 
sampling  problem  for  rainfall,  i .  e . ,  to  the  question  of  how 
many  raingages  are  required  to  achieve  a  certain  level  of 
confidence  in  runoff  prediction.  The  second  aspect  of  the 
spatial  pattern  problem  is  the  question  of  whether,  or  to 
how  great  an  extent,  spatial  variability  in  precipitation, 
in  and  of  itself,  affects  measurable  attributes  of  stream 
hydrographs  given  that,  in  an  actual  hydrological  record, 
other  sources  of  error,  uncertainty  and  spatial  variation 
may  obscure  any  relationship. 

The  influence  of  rainfall  variability  on  simulated 
streamflow  has  been  approached  primarily  with  the  rainfall 
sampling  problem  in  mind.  Dawdy  and  Bergman  (1969) 
concluded  that  errors  in  estimating  precipitation  volume  and 
intensity  over  a  catchment  were  likely  to  be  the  factor 
limitinq  the  accuracy  of  runoff  simulation  in  most  cases. 
°hanartzis  (1972)  demonstrated  the  importance  of 
incorporating  an  altitudinal  pattern  in  precipitation  in 
simulating  runoff  from  a  watershed  in  the  San  Dimas 
Experimental  Forest.  Wilson  et  a  1  .  (1978  )  compared 


simulated  hydroqraphs  generated  using  a  distributed  rainfall 
input  with  those  generated  using  rainfall  for  a  single  gage 
for  the  same  storm.  They  concluded  that  errors  in  the 
estimation  of  precipitation  input  may  result  in  serious 
errors  in  predicted  runoff  hydrographs.  On  the  other  hand, 
Chow  (1978)  concluded  that  a  single  rainqage  input  was 
sufficient  for  establishing  a  relationship  between  monthly 
precipitation  and  monthly  streamflow  for  the  Sanqamon  River 
in  Illinois. 

In  this  paper  we  examine  the  importance  of  spatial 
variability  in  precipitation  from  both  the  standpoint  of  the 
sampling  problem  and  the  problem  of  isolating  effects  of 
spatial  precipitation  patterns  frcm  a  noisy  record.  The 
problem  was  investigated  for  a  particular  case  study  -  a 
catchment  in  central  Illinois  -  as  part  of  the  V I N  project, 
a  meteorological  experiment  carried  out  in  central  Illinois 
by  the  University  of  Virginia,  the  JJlinois  State  Water 
Survey  (ISWS)  and  the  Rational  Hurricane  and  Experimental 
Meteorology  Laboratory.  Data  from  this  experiment  and  from 
earlier  experiments  conducted  by  ISWS  were  used  to  examine 
the  effects  of  rainfall  pattern  of  summer  storms  on  stream 
hydroqraphs.  A,  deterministic/statistical  modelling 
procedure  was  used  to  estimate  the  relationship  between 
length  of  record  available  and  the  level  of  confidence 
attached  to  the  differentiation  of  differences  in  hydrograoh 
characteristics  resulting  from  storms  of  different  patterns. 
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AVAILABLE  DATA 

Data  from  the  VI N  experiment  were  collected  during  July 
and  August  of  1979.  Rainfall  measurements  at  a  frequency  of 
5  minutes  from  an  array  of  gaqes  (see  F i q .  1)  spaced  at 
approximately  5  km  intervals  were  available.  Bihourly 
stream  flow  data  for  Friends  Creek,  a  catchment  of  287  km*-, 
were  provided  by  the  USGS.  The  Friends  Creek  catchment  is 
outlined  on  Fig.  1. 

Precipitation  data  for  the  East  Central  Illinois 
Network  (EC  IN),  which  was  operated  by  the  ISWS  durinq  summer 
periods  for  12  years  beginning  in  1955,  were  made  available 
to  us  by  ISWS.  Information  on  the  characteristics  of  the 
network  and  on  the  t i me/space  patterns  of  rainfall  observed 
in  the  EC  IN  network  are  given  in  Huff  (1967,  1968).  Stream 

gage  records  for  the  1955-1958  period  for  Goose  Creek,  a 
catchment  of  122  km^  within  ECIM,  were  obtained  from  the 
USGS.  \ 

THRESHOLD  ANAL YS I S  OF  EC  IN/GOOSE  CREEK  DATA 

The  threshold  analysis  used  in  this  study  is  an  attempt 
to  separate  the  effects  of  spatial  variability  of  convective 
precipitation  on  storm  runoff  from  the  many  other  influences 
on  catchment  response,  such  as  intensity  and  duration  of 
storm  omcioitation  and  the  pattern  of  antecedent  moisture 
conditions  over  the  catchment.  The  concept  is  to  separate 
the  events  in  the  rainfall  record  into  several  categories 


3 


based  on  the  pattern  of  storm  precipitation,  and  then  to  use 
si -pie  statistical  techniques  to  test  the  hypothesis  that 
rainfall  pattern  has  a  siqnificant  effect  on  the 
characteristics  of  the  resultant  hydrographs.  The  initial 
analysis  was  carried  out  using  the  ECI’l /  Goose  Creek  data  in 
the  summers  (June  to  August)  of  the  overlap  years  1355/53. 
'his  gave  a  sample  of  155  rainstorms  for  analysis. 

The  separation  procedure  used  is  shown  schematically  in 
Fiq,  2.  The  procedure  accepts  average  storm  precipitation 
from  three  subcatchments  of  the  basin.  On  the  basis  of  a 
catchment  a r ea /d i s t a n ce  from  outlet  histogram,  three 
subcatchments  of  equal  area  were  defined.  Average  storm 
precipitation  in  each  subcatchment  was  calculated  by 
multiplying  the  gage  totals  by  appropriate  Thiessen  polygon 
weighting  factors.  Suitable  values  of  the  threshold 
parameters  were  found  to  be  T H 1  =  1 . 2 5 ,  T  H  2  =  5 . 0  mm,  TH3  =  0.2 
and  TH 4  =  5 . 0  mm.  It  was  found  that  storms  showing  distinct 
spatial  pattern ns  in  total  storm  rainfall  (the  temporal 
distribution  at  each  site  was  not  available)  were  relatively 
rare.  Only  5  useable  storms  were  separated  into  category  1 
and  5  in  cateaory  2.  These  storms,  together  with  the 
characteristics  of  the  associated  hydrographs  ere  shown  in 
Table  1.  Given  this  sample  size  none  of  the  characteristics 
were  statistically  different.  Mean  values  of  a  v  e  r  a  a  e 
catchment  precipitation,  delay  time  and  rise  time  for  the 
two  categories  were  very  close.  Any  differences  that  there 
•’v  be  in  "■?  response  of  the  two  categories  are  obscured  by 
differences  in  other  factors. 
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TAlll.E  1.  Soloclcd  rainfall  and  hydrograph  charac  ter  in  Lies  for  tort.ir;  in  wh 
the  precipitation  v/as  highest  on  the  upper  basin  (Category  1)  and 
those  for  which  the  highest  precipitation  was  on  the  lower  basin 
(Category  2). 


o  ■— 
n 

•H  O 
ei  2 


in 

n  (N  o  rn  (N 


< 


O  m  O  CN’  rj 

•  •  •  •  • 

n  o  o  N  n 


S> 

Si 


CO 

CA 

CN 

r- 

o 

c 

c 

rH 

CN 

CA 

r- 

a 

c~> 

• 

• 

• 

• 

• 

cs 

G 

CN 

CN 

m 

r—{ 

m 

si 

10 

s> 

s> 

c 

o 

— ' 

GJ 

NO 

o 

O 

CN 

a 

a 

• — i 

n 

in 

o 

o 

in 

• 

# 

• 

• 

• 

o 

G 

»H  i — { 

CN 

• — i 

CN 

si 

Si 

c 

co 

< 


0  r-I 

Cn  G 
G  si 
O  C 
O 

>  G 

<  a 


rl  O  Mn  C3 

in  n  o  h  o 


o  o  O  r-l  o 


i-i 

0 

0  o 

■o  0 
g  u 
c  n 


r->. 

in 

o 

r- 

p-. 

CO 

>1 

in 

m 

m 

m 

LO 

o 

0 

r- 

c 

/-■ 

r— f 

? 

cr 

3 

o 

n 

< 

s 

c 

m 

CN 

CO 

n 

U 

CN 

r— i 

r-i 

in 


■r  -r  H  r:  n 


on i  o  cn  o 
o  --I  o  'r 


<o  in  ci  o  ro 
ci  *=■  o  r~ 


<o  o  m  o  co 


Mn  n  o  <? 

in  CO  "C  C3  CM 


OONNCO 


ro  rs  "C  in  ni 

IT?  CN  f'  CN 


o  o  o  o'  o 


<N 

o 

o 

CO 

CO 

& 

in 

in 

LO 

in 

LO 

0 

cr. 

£ 

£ 

u  1 

o 

o 

3 

o 

< 

< 

‘*2 

•  o 

m 

o 

O 

rj 

n 

*■— ? 

r-i 

* — { 

Si 


0 

o 

n 

a 

•H 

2 

S-l 

•H 

S 

<- 

c* 

O 

G 

Cl 

tl 

G 

Cl 

U 

0 

G  • 

!-( 

>  s> 

'w 

G  G 

>1 

O 

o  s 

a 

•o  o 

-U 

si 

G 

n  g 

M 

3  o 

r* 

O 

■H  O 

— 

E  -C 

s3 

SI 

CJ 

r-r 

Cl  c 

•H 

G  -h 

S)  JC 

o 

11  SI 

in 

• 

•H 

2 

G 

C  > 

G 

tn 

•  f-* 

G 

11 

O 

s) 

O  O 

w 

n  o 

o 

ns  g 

r— 1 

G 

o  t7> 

n 

•H 

n  c 

G 

U  H 

O 

G 

G  G 

■H 

11 

•H  W 

SI 

G 

O 

G 

i— i  si 

.r— 

'r-l 

G  G 

*  *""" 

O 

•H 

S) 

G 

s)  G 

n 

•H 

■H  O 

u 

G  -iH 

•H 

r — 1 

••H  S 

si 

G 

G 

•H 

SI  S 

S| 

Si 

O  -H 

O 

•H 

a. 

C 

o  -H 

O 

•  H 

5  o 

E 

•h  a 

•H 

S| 

s)  o 

O 

0, 

o 

O 

G 

s:  >— i 

r; 

2 

Sl  o 

U 

•H 

SI 

n  c 

W 

G  O 

G 

•H 

■— 

n-  si 

_> 

a  g 

a 

G  H 

i.1 

■r-i  S 

•r-i 

3 

si  -H 

1 

O  G 

CJ 

•S 

a  -h 

G 

2 

< 

< 

* 


7 


A  SIMPLE  DISTRIBUTED  RAI NF RLL/RUMQFF  MODEL 


The  discussion  above  suggested  that  given  the  data  base 
available  from  the  ECIN/Goose  Creek  measurements,  it  would 
not  be  possible  to  draw  any  firm  conclusions  regarding  the 
importance  of  Knowledge  of  the  spatial  variability  of 
rainfall  inputs  to  predicting  storm  discharges.  The 
implication  is  that  the  four  summers  of  ra i n fa  1 1 /d i s cha rqe 
data  available  for  the  Goose  Creek  catchment  are  not 
sufficient  to  distinguish  the  effect  of  rainfall  pattern 
from  other  factors  such  as  antecedant  moisture  conditions, 
storm  durations  and  intensities.  Consequently  a  somewhat 
modified  question  was  examined.  How  long  a  period  of  record 
would  be  necessary  before  the  effect  of  spatial  variability 
could  be  demonstrated  to  be  statistically  significant?  To 
make  a  first  estimate  of  the  length  of  data  required  a 
deterministic/statistical  simulation  approach  was  adopted. 

The  steps  in  the  analysis  are  shown  in  Fig.  3.  The 
validity  of  the  analysis  depends  largely  on  the  quality  and 
quantity  of  the  data  available  for  calibration  of  the 
rainfall  and  rainfall/runoff  models  used.  The  best  data 
available  for  the  rainfall/runoff  calibration  were  the 
VI N /Friends  Creek  data  collected  durinq  the  summer  of  1979. 
Five  minute  rainfalls  and  bihourlv  streamflow  data  were 
available.  uowever,  the  discharge  record  durinq  this 
relatively  drv  summer  was  dominated  by  a  single  storm 
hydro graph,  preceded  by  two  smaller  hydrographs  as  the 
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catch  nent  wetted  up  durinq  the  end  of  July  and  beginning  of 
August.  The  previous  hydrograph  peak  of  any  significance 
was  on  April  26,  1979. 

It  was  clear  that  this  data  base  would  not  support 
calibration  of  any  complex  rainfall/runoff  model.  A  number 
of  simple  deterministic  distributed  models  were  investigated 
but  the  final  choice  was  the  simplest  of  all.  This  model  is 
shown  schematically  in  Fig.  A  and  consisted  of  a  number  of 
subcatchment  models  contributing  to  a  constant  wave  speed 
(linear  channel)  channel  routing  component  similar  to  that 
used  by  Surkan  (1969)  and  K  i r  k  by  (1976). 

The  subcatchment  model  consists  of  two  storage 
elements.  The  store  SI  is  associated  with  two  threshold 
parameters:  one  (Sc)  controlling  discharge  from  a  riparian 

contributing  area  Ac;  the  other  (Sd)  controlling  input  to 
the  S2  store.  Various  non-linear  forms  were  tried  for  the 
S2  store,  but  early  runs  showed  that  a  simple  linear 
formulation  was  satisfactory.  Evaporation  was  taken  from 
the  SI  store  at  the  measured  pan  potential  rate  (available 
on  a  daily  basis)  until  the  store  was  empty. 

The  channel  routing  procedure  is  equivalent  to  a 
constant  time  lag  for  each  incremental  length  of  the  channel 
network.  Input  per  unit  area  in  each  subcatchment  is 
assumed  constant  at  each  time  step.  Total  catchment 
discharge  is  then  given  by  the  discrete  linear  convolution 

*■  M 

V  V 

0  ( t )  =  Zj  Zj  d  ( t  -  i  .  j  )  a  ( t  -  i  ,  ,i  ) 
i  =  0  .1  =  1 
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where  Q(t)  is  catchment  discharge  at  tine  t,  M  is  the  number 
of  subcatchments,  q  ( t  ,  j  )  is  the  predicted  discharge  from 
subcatchnent  .i  at  time  t,  and  a  ( t ,  j )  is  the  fractional 
catchment  area  in  subcatchnent  j  having  a  travel  time  t  to 
the  catchment  outlet. 

The  number  of  subcatch  rents  is  arbitrary  and  could  have 
been  made  equal  to  the  number  of  raingage  input  records 
available.  However,  three  subcatchments  have  been  used  in 
the  present  study  to  be  compatible  with  the  threshold 
procedure  described  above.  It  was  felt  that  three 
subcatchments  would  retain  sufficient  r ep r e s e nt a t  i  on  of  the 
spatial  pattern  of  precipitation  but  has  the  advantage  that 
the  computing  requirements  in  running  the  model  ’during  the 
calibration  phase  are  not  tco  great,  and  specification  of  a 
stochastic  rainfall  simulation  model  is  greatly  facilitated. 
The  three  subcatchments  were  chosen  to  be  equal  in  area  and 
bounded  by  lines  of  equal  distance  along  the  channel  network 
(Fig.  5).  The  definition  of  the  area  contributing  to  each 
channel  reach  was  very  subjective  in  such  flat  terrain,  but 
is  not  crucial  to  the  analysis. 

CAl  IBRATIQN  OF  THjE  R  A  I  NF  ALL /RUA'OF  F  MODEL 

A  period  of  450  hours  was  used  for  model  calibration, 
’he  original  rainfall  records  were  lumped  into  hourly 
subcatchment  averages  usino  Thiessen  polvoon  weiohts.  A 
model  time  step  of  one  hour  was  used.  The  model  has  ‘our 
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parameters  to  be  estimated  for  each  subcatchment  and  one 


channel  routing  parameter.  In  addition  the  initial 
conditions  at  the  start  of  the  simulation  for  each 
subcatchment  must  be  specified.  For  simplicity,  the 
parameter  values  were  assumed  to  be  equal  for  all  three 
subcatchments.  After  the  long  dry  period  preceding  the 
calibration  period  the  store  SI  was  assumed  to  be  empty. 
Outflow  per  unit  area  from  the  S2  store  was  set  equal  to  the 
observed  catchment  discharge  at  the  first  time  step  in  the 
record  for  all  three  subcatchments.  In  addition,  early 
trials  suggested  that  the  riparian  contributing  area 
parameter  Ac  could  be  fixed  at  its  initial  estimate  of  0.02 
(the  estimated  fractional  area  of  channel  bed  and  banks). 
This  left  only  four  parameters  to  be  calibrated  (Sc,  Sd,  K 
and  Cv)  and  minimized  parameter  interaction. 

The  parameters  were  calibrated  using  a  Rosenbrock 
automatic  optimization  procedure  available  from  the 
Institute  of  Hydrology  modelling  packaqe  (Oouqlas,  1974). 
The  optimized  values  are  given  in  Table  2  and  the  observed 
and  simulated  hydrographs  for  the  calibration  period  are 
shown  in  Fig.  5.  The  final  objective  function  value  gave  a 
modelling  efficiency  F  of  96.4%  where 


450 

E 


F  -  !  _  111 

no  "  1  450 


(  Q05S  ( t )  -  QSIH.(t))2 


(Q0BS(t)  -  Q) 


t  =  i 


and  Q  is  the  mean  observed  discharge,  9  0  a.  S  ( t )  is  trie 
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observed  discharge  and  QSIM(t)  the  simulated  discharge  at 
tine  t.  Given  the  lenqth  of  record  available  there  was  no 
possibility  of  doinq  a  split  record  model  validation  in  this 
case. 

The  optimized  channel  wave  velocity  in  Table  2  may  be 
compared  with  measured  tracer  travel  times  (givinq  f 1 o w 
velocities)  in  the  upper  23  km  of  the  Sangamon  River  (of 
which  Friends  Creek  is  a  tributary)  in  Table  3.  A  channel 
wave  velocity  of  0.597  km/hr  gives  a  maximum  wave  travel 
time  from  the  furthest  headwater  of  41  hours  and  a  mean  wave 
travel  time  (defined  by  the  centroid  of  the  catchment 
area/travel  tine  distribution)  of  19  hours. 

The  fact  that  the  spatial  pattern  of  precipitation  may 
be  important  in  the  storm  runoff  process  is  suggested  by 
running  the  model  using  the  optimized  parameter  values  but 
with  spatially  homogeneous  inputs.  Fig.  7  compares  the 
observed  discharge  hydrograph  with  those  predicted  using  the 
area  A  input  all  over  the  catchment  and  area  C  input  all 
over  the  catchment.  These  predictions  are  markedly 
different  from  the  spatially  distributed  input  prediction  of 
Fig.  6.  Of  course,  the  model  might  have  been  calibrated 
using  only,  say,  area  A  input  and  a  reasonable  fit  obtained 
given  the  relatively  short  lenqth  of  record.  Our  results 
still  demonstrate  that  the  effect  of  using  a  non- 
representative  set  of  gages  (relative  to  calibration)  is  of 
considerable  importance. 
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Calibrated  values  of  model  parameters. 


CATCHMENT  ? A RA.METERS 


Fractional  riparian  contributing  area  0.02 


SI 

store  threshold 

for  riparian  area 
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for  input  to  S2 

105. 

.  60  mm 

me 

an  residence  tin 

.e  in  S2  store 

56, 

,  8  2  hr 

.'.CP EL  ROUTING  PARAMETER 

Wave  velocity  0.5S7  km/hr 


3.  Measured  channel  flow  velocities  (after  Stall 

and  Hiestand,  1969,  for  upper  23  km  of  Sangamon 
River,  origin  to  Saybrook) . 
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RAINFALL  MODEL 


The  11-vear  rainfall  record  from  EC  IN  was  used  to 
develop  a  procedure  to  generate  a  synthetic  record  of 
precipitation  events  with  statistical  properties  appropriate 
for  central  Illinois.  The  E  C  IN  network  was  superimposed 
over  the  F-iends  Creek  catchment  and  Thiessen  weights  were 
used  to  compute  rainfall  on  the  three  subareas  A,  8  and  C. 
Average  depths  and  durations  for  each  of  the  three  areas  in 
all  summers  (June,  July  and  August)  were  determined.  Storms 
with  more  than  3  showers  or  with  duration  greater  than  24 
hours  or  with  depths  less  than  0.1"  were  eliminated  fron  the 
analysis. 

A  number  cf  stochastic  models  for  thunderstorm 
synthesis  have  been  developed  (e.g.,  Sorman  and  Wallace 
1972;  Duck  stein  et  al.  1972  ;  Smith  and  Schreiber  1974  ; 

Croley  et  al.  1978).  Rather  than  employing  a  sophisticated 
stochastic  rainfall  generation  model,  we  chose  a  simple 
statistical  procedure  as  consistent  with  the  goals  of  our 
study.  Storm  events  were  generated  independently  and  the 
statistical  effect  of  different  interstorm  time  periods  was 
accounted  for  by  choosing  initial  conditions  for  the  S2 
store  of  the  runoff  model  from  an  exponential  distribution 
represent! no  the  range  of  initial  (pre-storm)  discharges 
observed  o  the  Goose  Creek  and  Friends  Creek  records 
-  ■;  s  u  r  e  9 ) .  Leckim  further  information  the  initial  value 
of  the  SI  store  was  taken  from  a  uniform  distribution  with 
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Figure  8.  Frequency  histogram  of  initial  discharges  per  unit  prior 
to  storm  hydrographs.  Combined  Goose  Creek  and  Friends 
Creek  (shaded)  data. 
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the  r  a  n  q  e  (  0  ,  S  „• ) . 

The  rainfall  generation  scheme  must  produce  1)  total 
deocn,  2)  duration,  3)  spatial  pattern  and  4)  temporal 
distribution.  The  first  three  items  were  selected  from  a 
multiv  ariate  log-normal  distribution.  A  scatter  diagram  or 
depths  and  durations  of  storms  for  area  3  suggested  that  a 
logarithmic  transformation  of  the  data  would  be  appropriate. 
Fic.  9  shows  the  frequencies  of  log  values  of  these 
quantities.  Spatial  pattern  was  preserved  by  the  use  of 
depth  "multipliers'1  for  A/3  and  C/3.  That  is,  the  ratios  of 
storm  depths  for  area  A  to  area  B  and  area  C  to  area  3  were 
calculated  ‘‘or  each  summer  storm.  Rainstorms,  including 
spatial  pattern,  can  then  be  simulated  (statistically)  by 
choosing  depth  and  duration  for  area  3  and  the  A/3  and  C/B 
multipliers  from  an  appropriate  distribution.  Depth  for 
area  A,  for  example,  is  then  calculated  simply  as  the 
product  of  multiplier  A  and  the  depth  for  area  B.  A  sample 
of  the  observed  storm  multipliers  is  plotted  in  Fig.  10  and 
the  histograms  of  log  values  are  shown  in  Fig.  11.  The 
means  and  the  variance-covariance  matrix  for  log  values  of 
depth,  duration,  multiplier  A  and  multiplier  C  are  listed  in 
Table  4,  Although  the  hypothesis  that  the  storm  depth  and 
the  multipliers  came  from  a  loq-normal  distribution,  can  be 
rejected  at  the  95%  level  by  a  Chi  square  goodness-of-fit 
test,  tract i cal  considerations  for  generation  of  correlated 
r a ' : c ~  variables  dictated  use  o'  a  multivariate  lea-normal 
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MULTIPLIER,  AREA  C  (LOG) 

Figure  11.  Frequency  histograms  of  log  values  of  storm  multipliers 


distribution.  e  do  not  think  that  our  general  conclusions 
depend  strongly  on  this  assumption. 

Tine  distributions  of  rainfall  were  derived  from  the 
normalized  cumulative  precipitation  curves  for  thunderstorms 
presented  by  Huff  (1967).  Storms  were  assigned  to  one  of 
the  four  types  of  hyetographs  shown  in  Fig.  12  according  to 
the  following  probabilities:  type  I,  0.37;  type  II,  0.39; 
type  III,  0.20;  type  IV,  0.04. 

Daily  evaporation  was  modelled  as  a  normally 
distributed  random  variable  independent  of  rainfall.  Data 
for  summer,  1979  showed  no  significant  difference  between 
evaporation  on  rain  v_s_  no-rain  days. 

SIMULATION  RESULTS 

One  thousand  storms  were  generated  using  the 
statistical  procedure  described  above.  This  represents  a 
sample  of  approximately  75  to  100  years  based  upon  the 
observed  E  C I N  network  statistics.  The  threshold  analysis  as 
previously  described  was  applied  to  the  synthetic  record  and 
storms  were  classified  as  either  concentrated  on  the  upper 
basin  (category  I),  concentrated  on  the  lower  basin 
(category  II)  or  relatively  uniform  (category  III). 
Hydrographs  were  simulated  for  the  storms  in  each  class  and 
peak  discharge,  time  to  peak  and  total  storm  volume  were 
noted. 

The  potential  for  isolating  the  effects  of  spatial 
rainfall  patterns  from  other  : nee r t a i nt i es  can  be  gaged  by 


TABLE  4.  Statistics 

for  leg  values  of 

depth,  curat 

ion , 

multiplier 

A  and  multiplier 

C. 

■lean  Values 

Depth 

-1 .084 

Duration 

1.104 

V.-q  T  z~  **  D*)  '  ar 

-0.093 

Multiplier  C 

-0.144 

Var iar.ce-Ccvariance  .’ 

latrix 

Depth  Duration 

Multiplier  A 

Multiplie 

Depth 

0.637  0.355 

0.029 

0.057 

Duration 

0.810 

0.109 

0. 061 

Multiplier  A 

0.276 

-0. G99 

Multiplier  C 

0.326 

TABLE  5.  Kolmogorov-Smirnov  statistics  comparing  hydrograph 
characteristics  for  the  three  storm  categories 
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m ,  n 


Category 

Category 

Category 

I  vs  1 1 

II  vs  III 

I  vs  III 

Peak  : lev 

0 .136 

0. 17S 

0.267 

Time  to  Peak 

0.767 

0.294 

0. 506 

Total  Stcrmficv;  Volume 

0.073 

0.037 

0.118 

9  5 %  value  of  d 

0.176 

0.193 

0.199 

to  vo  cr 

Ov  VD  * 


three  threshold 


comparing  hydrograph  characteristics  for  she 
categories.  Table  5  gives  the  values  of  d  n  n,  the  two- 
sample  Kolmoqorov-Smirnov  statistic,  for  the  three 
categories  and  for  the  three  hydroqraph  characteristics. 

The  value  of  d  required  for  refection  of  the  null  hypothesis 
at  the  95 %  level  is  also  noted.  Time  to  peak  discharge  is 
significantly  different  amonqst  cateoories  but  except  for 
category  I  peak  discharge  vjs_  category  III  peaks,  no  other 
differences  are  significant  even  with  ca.  100  years  record. 
Given  that  records  of  length  tens  of  years,  let  alone  100 
years,  are  not  readily  available  for  the  type  of  problem  we 
are  consider'nq,  a  precise  answer  to  the  question  of  how 
long  a  record  is  required  was  not  sought,  it  suffices  to  say 
that  a  very  long  record  indeed  would  be  necessary. 

To  compare  the  importance  of  spatial  rainfall  patterns 
per  se  with  the  importance  of  the  samplinq  problem 
associated  with  rainfall  variability,  hydrograph  simulations 
were  also  prepared  usinq  area  A  and  area  C  precipitation  as 
uniform  inputs  to  the  runoff  model.  That  is,  hydrographs 
were  simulated  as  if  only  data  from  area  A  (or  area  C)  were 
available  and  as  if  those  data  were  extrapolated  to  the 
entire  catchment.  Table  6  lists  values  for  the  Kolmogorov- 
Snirnov  statistic  for  comparisons  between  distributed  input 
simulations  and  those  done  with  the  input  approximated  as 
uniform  based  on  either  area  A  or  area  C.  Differences 
amongst  peak  discharoes,  time  to  peak  and  total  v  o  1  u  -  e  are 
a  1  ’  significantly  different  'or  c a t e n o r v  !  ,  distributed  v_s_ 

area  A  input,  and  for  cateqory  II,  distributed  vs  area  C 
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TABLE  6.  Kolnogorov-Smirnov  statistics  comparing  hydrograoh 
characteristics  for  different  inputs  in  various 
categories . 


d 

m ,  n 

Spatially  distributed  Spatially  dis 

input  v_s  input  for  input  vs  inpu 

area  A  area  C 


CATEGORY  I 


Peak  flow 

0. 

Tima  to  peak 

0. 

Total  volume 

0. 

9  5%  d 

0. 

m,n 

9  C  £ 

0.139 

74  9 

0.472 

261 

0.227 

197 

0.182 

CATEGORY  II 


Peak  f lev/  0.094 

Time  to  peak  0.445 

Total  volume  0.143 

95%  d  0.170 

m,n 


0.24  5 
0.349 
0.243 
0.182 


CATEGORY  III 

Peak  flow  0.074 
Time  to  peak  0.111 
Total  volume  0.086 
95%  d  0.214 


0.071 

0.079 

0.070 

0.215 


-IJ  4> 


incut.  As  expected,  none  of  the  differences  were 
siqnificant  for  the  cateqory  III  storm  hvdroqraphs.  These 
results  are  for  ‘he  entire  sequence  of  storms. 

A  similar  ccmoarison  between  cateqory  I  storms  - 
distributed  vs_  area  A  input  -  and  catecory  II  storms  - 
distributed  v_s_  area  C  input  with  1/2  of  the  record  length  is 
given  in  Table  7.  Again  all  characteristics  show  highly 
significant  differences  and  we  infer  that  even  smaller 
samples  would  also  show  siqnificant  differences. 

DISCUSSION 

The  results  from  the  analysis  of  the 
deterministic/statistical  modeling  data  indicate  that 
differences  in  the  timing  of  peak  flows  are  most  likely  to 
be  discerned  as  an  effect  of  spatial  rainfall  pattern. 
Differences  in  peak  discharges  are  much  less  likely  to  be 
isolated  from  a  background  of  other  influences. 

Un f or t u na t e 1 y ,  in  real,  as  opposed  to  synthetic,  streanflow 
records,  Peak  flows  and  their  timing  can  be  difficult  to 
define  unambiguously  because  of  complex 

hyetograoh/hydrograph  shape.  This  fact  may  make  the  problem 
of  discerning  effects  of  rainfall  pattern  even  more 
difficult  than  suggested  bv  this  study. 

On  the  other  hand,  storm  movement,  which  was  neglected 
in  ,  r  s  *.  j  d  v ,  might  accentuate  effects  of  spatial  pattern  in 
s  o " e  "nstances.  for  the  -ri e-  is  Creek  catchment  differences 
in  *ne  t  ninq  of  summer  thunderstorms  are  small  compared  to  , 


Kol~.ogorov-Snirn.ov  statistics  comparing  hydrograph 
characteristics  for  category  I  storms  (distributed 
vs  area  A  inputs)  and  for  category  II  storms 
(distributed  vs_  area  C  inputs)  for  half  of  the  total 
record  length. 


TZGORY  I 

Peak  flow 
Time  to  peak 
Total  volume 


Distributed  vs  area  incuts 


0.316 

0.772 

0.305 

0.279 


.TZGORY  II 

Peak  flow 
Time  to  peak 
Total  volume 


0.336 

0.333 

0.276 

0.257 


31 


the  channel  tine  of  travel  and  consecuently  our  results  are 
unlikely  to  be  affected  by  sophistication  of  the  rainfall 
modellinq  procedure.  We  have  also  neglected  several  other 
factors  that  night  accentuate  the  effects  of  rainfall 
pattern  under  certain  conditions,  for  example  spatial 
differences  in  catchnent  response,  spatial  differences  in 
antecedent  conditions  and  diffusion  within  the  channel 
network.  Much  nore  comprehensive  discharge  records  and 
channel  data  w o u 1 d  be  required  before  such  effects  could  be 
studied.  Our  results  suaqest  that  only  in  cases  where  such 
effects  were  particularly  marked  would  rainfall  pattern  have 
a  significant  effect  on  peak  flows  and  storm  runoff  volumes. 

Although  the  effects  of  rainfall  pattern  per  s e  may  not 
be  striking,  the  analysis  of  the  sampling  problem  indicates 
that  significant  differences  in  predicted  peak  discharge 
will  result  in  cases  where  a  "non-representative"  subsample 
of  raingages  is  used  to  provide  input  to  the  model.  Thus, 
we  conclude,  as  did  Dawdy  and  Bergnan  (1959)  and  Wilson  et 
a  1  .  (1978),  that  an  accurate  portrayal  of  spatial  variation 
in  rainfall  is  a  prerequisite  for  accurate  simulation  of 
streamflows. 
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